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ABSTRACT 



The diabatic heating rate at lOOmb is determined from the thermo- 
dynamic equation, utilizing previous calculations of the potential 
function x and the n vertical velocity” w, for the period 1-5 April 1963. 
The four basic terms of the expanded thermodynamic equation; storage, 
stream-velocity, x~P°tential advection, and vertical heat transport are 
examined individually, and added to obtain Q, the diabatic heating rate 
which is compared to the results of previous investigations. The 
zonally-averaged heating rates are discussed in terms of stratosphere 
circulations and their associated physical mechanisms. 
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1 . 



Introduction 



There has been recent interest in the diabatic heating rates and 
related energy changes in the upper troposphere and lower stratosphere, 
notably by Davis [1], Julian [2], Martin [4], K. Myakoda [3], and others, 
Diabatic heating is, of course, also at the heart of the problem of 
large-scale energetics. However, the usual approach through energetics 
yields average values for a hemisphere for the period of study, whereas 
the approach of this study takes the more direct route of obtaining 
12-hour average heating rates at each point of the usual octagonal 
(51 by 47) analysis-prediction grid system employed by the Joint Numer- 
ical Weather Prediction Unit, and other groups forecasting large-scale 
weather changes. 

Martin [4] demonstrated the feasibility of the grid-point approach 
by obtaining 12-hour mean values of both the potential function^ of 
the wind and the "vertical velocity" Once these fields are 

known it is possible, utilizing the thermodynamic equation, to solve 
for Q, the diabatic heating rate. Martin’s model is described in his 
paper [4] and only those aspects required for this paper will be re- 
viewed here. For more details regarding Martin’s model, the reader is 
referred to his paper. 

The data period selected for the calculations was 1-5 April 1963. 

The data was provided by the Joint Numerical Weather Prediction Unit, 
Washington, D.C. This particular data period was selected in an effort 
to reflect the relatively large heating rates and the associated large 
meridional flows that normally occur during April, particularly in the 
polar stratosphere. 
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2. The diabatic equation 



If one expands the equation of thermodynamics, the result 




follows. Here Q is the heating rate per gram, while ^ and X give 
the nondivergent and irrotational parts of the wind in accordance with: 

V^fexv^+vX 

Here the fields of ^ and X were obtained preliminarily by Martin [4] 
by solution of the Balance Equation^ and the x~2 equation respectively 
(see section 3). When the equation of state and hydrostatic equation 
are employed in (1), we obtain our basic working equation: 




The fourth term of (2) may be transformed into the more convenient form: 



& / 

/ (\oo) 

ap 



n\i% 




(3) 



where <T = - but is also the negative of the coefficient of a ) 

9 v , V 

in the final form of (3). Here { and ( are centered at 75mb 

\^p/^ Up# 

and 125mb respectively, that is, at the upper-half and lower-half layer 
centers, respectively. The term - (T* of (3) has a positive sign since 
both /<^p and the contents of the last bracket are negative. This 
agrees with the intuitive expectation that an updraft (w<o) is associated 
with vertical-motion cooling. 

The factor is present in all terms on the left side of (2), 
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and in order to determine its value at lOOmb, the 50mb and 150mb data 



levels must be utilized. Unfortunately there was no vertically con- 
sistent 50mb data in punched-card form immediately available, so it was 
necessary to derive the 50mb contour data completely from regression 
equations. In order to do this the United States Navy Weather Research 
Facility Regression equations of D. Lea [5] were used. These equations 
have the following form: 



where K Q , K^, K£ represent the regression coefficients. It should be 
noted that the regression coefficients K Q , K^, K 2 are listed by ten 
degree latitude bands from ten degrees north to the pole. It was 
necessary, therefore, to devise a computer program converting each 
K n (‘| } ) to a value of ^ = K n ( i , j ) defined at each grid point within 
the 51 by 47 octagon. This conversion was based upon a subprogram 
giving sin (j) in terms of the (i,j) grid. In addition K n (^) was made 
continuous with respect to latitude by averaging the values across the 
defining latitude boundaries for each ten degree band, while retaining 
the latitude^band mean value. 

Since the lOOmb contour analyses during the period of study were 
considered to be geographically consistent, whereas, those of tempera- 
ture were not, it was decided again to obtain T^qq using T-^qq = -A 



in (4). With this modification, the regression equation is expressible 
in the following form: 






( 5 ) 
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where 



i 



K 

K 



/ 

1 

/ 

2 



(g 50 “ 2 100 ) " 

and the term 



^ ^ p 



(Kj_ - 1) 




estimated thickness of the layer 100 to 50mb. 
appearing in ( 2 ) has been expressed in the form: 



— C^5o *-ioq) fcioo 

JU3 




( 6 ) 



In (5) and ( 6 ), (#50 - 3^qq) is estimated from the modified regression 
equation (5), with input data 2 ;loO anc * (**100 ” ^150^ obtained from 
copies of data tapes provided by the Joint Numerical Weather Prediction 
Unit. The parameters £;lq 0 end - ^ 59 ) of (5) and ( 6 ) were used 

wherever possible in order to utilize as much vertically consistent 
data as possible, thereby tending to minimize errors due to the use of 
the regression equations (4) or (5). All computations were performed 
using the CDC 1604 computer. 



Equation (3) can now be written in its final form: 




It should be noted that 3 = D + where «frp is the standard atmosphere 
height and D = £■ - is the height anomaly relative to that of the 
standard atmosphere. In the first three terms of (7) 2p is independent 
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of the variables x,y,t so that the time and space derivatives become 
simply: 




etc., and hence 2*p can be omitted. However this cannot be done in the 
vertical motion term where affects the vertical stability. 

3. The irrotational wind and vertical motion structure. 

As pointed out in section 1, the development in this section follows 
Martin [4], who solves diagnostically for X for his X“2 equation which 
is written: 




This equation is solved by relaxation using an iterative, scan-procedure, 



subject to the convergence requirement: 

I *4 ~k\ 



< 4.0 



cm 



(9) 



'Lj 'Hj 

Then with $ assumed to be known, one may compute 

v- V or -|V^( = - (ty) 

Using m = G o (p) , as a cubic polynomial in p, with an upper boundary 
condition oj = 0 at p = 0, Martin also derives an integrated form of (9) 
giving a) = at each grid point. 

In terms of finite difference notation w (100) has been shown [4] 



to be: 



(ido) 



10 ; 



l j 




( 10 ) 



Here x is scaled in terms of centimeters^ as this scaling was defined 



5 



is ( 8 ) . From (10) , oj 



( 100 ) 



ij 



(lOo) 



takes on the finite difference form: 
2 



uo, 



_ cy rm_ 

~ ^ A 7 - 
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( 11 ) 



with assumed zero on the two outermost rows of the octagon. This 

was also the boundary condition used by Martin in solving ( 8 ) . The 

47 by 51 octagonal grid mapping program used in the computations was furnished 

on tapes by the National Center for Atmospheric Research. Since it is 

not possible to compute^^^f on the "outer boundary" (Fig. 1), it was 

necessary to set co = o on the outer boundary throughout the computations. 





The normal five-point map grid 1-2-3-4-0 was used one row within the 
computational boundary, while a nine-point grid was used at all points 
further inside. The nine-point grid simply inserts a three-point mean 
for a normal difference; for example, the X 3 “ Xi difference of a five- 
point grid is replaced with the weighted difference: 

in the nine-point grid. The actual five-point finite-differencing 
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schemes used in the computations are as follows. 




relcltt 



to tdf rrj 



\ r. 





Furthermore the resultant from (11) have been smoothed twice every- 
where interior to the two boundary rows just mentioned, giving , the 



values used in this study, 
4. Results 



With the fields of Xij and known, as well as and observed 

values of Ci A 1 la \Kr f ~rt \ 

M5 



I (z 



from consecutive 12-hour map-fields, Q has been determined from (7) 
The first term of (7), called the f, storage term” is 









l 



( 12 ) 



The second and third terms represents the 12-hour mean horizontal heat- 
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transport at the gridpoint (i,j) and have the combined form: 




( 13 ) 



The fourth term in (7) is the vertical heat transport (u>-transpor t) 




(14) 



Th6 four terms on the left side of (7), which have been listed in 
the preceding paragraph, have also been computed individually as 
functions of (i,j) everywhere within the computational boundary. Com- 
bined values have been formed to give resultant octagonal fields of 

• # 

Q(i,j). For a typical field of Q see Fig. 3. Zonal-averages of each 
of the five fields have been computed by programming an averaging 
routine for each five degree latitude band (over an area of 20-90 
degrees latitude) . 

Q, and the preceding four terms as well, have units of gpm which, 
for later comparison purposes must be changed to langleys. Q may be 
determined by a sum of terms of the form: 



^ "" ^ 4 100 



M c. p 3 cL(A2) 
Rd (U3^t 100 



(15) 



♦ — 1 1—9 

Here Q is in ergs(gm) 1 (sec)" 1 cm z provided all terms on the right 
side of (15) are in c.g.s. units , in gp-cm and M-^qq = mass of 
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90E 




Fig. 3. Example of printout of Q values for (06 GMT, 1 April 1963) 
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Oo 




the column 150mb to 50mb. 



(16) 



The expression within the braces of (16) is now identical to (4T/4t) 
in units of (°K)(sec) - ^. To get a corresponding value of Q per gram 
per day, both sides of (16) are multiplied by 86,400 sec (dayj} and 
divided by 102j the latter being the number of gms in lOOmb: 



q(JH —\ =. 

^ day) 



Sf. 

Wb 



3.41 1 jqcfx 8£40o lll-( A ^l 

IOZ Jb^ J 



= I34 1 1 q x 1 6 \ 84 4 ooIi(A £ ) 

i-cflS4L 1aTjl3t 



3?% 



sec 



3py : 



sec 



Q = 6.298 x 



[ % (4z -3 



ly (gm) 1 (day) 1 



The units of ^ /^t [A 2-] are emphasized as gpm (sec) ^ because this is 
the form of the computer output units. Thus it follows from 
that a geopotential change of 1 gpm (gm) “^sec”! is equivalent to a 
heating rate of: 



6.298 ly (gm) 1 (day) 



>-l 



The last statement is typical of each of the terms in (12), (13), (14). 

# 

Zonal averages for each term and their zona l mean sum [Q] of (7) have 
been determined and are listed by latitude belts in Tables 1 through 5. 
This is done also, for purposes of comparison, with computations of [Q] 
made by Davis [1], who used a climatological or standardized atmosphere 
(Table 6). Note that the notation [ ] indicates a zonal average at a 
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specified latitude. 

In Table 6, the nine sets of results obtained here are averaged 
to give a mean daily average of the storage, the horizontal transport 
terms of both the stream and potential types, the vertical-heat 

• _ i _ -I 

transport and the resultant Q, each in ly(gm) day. Davis 1 results 
have been centered at 117. 5mb, but presumably this slight shift of 
levels is insignificant. In regard to the units of ly(mb)“^ day”'*' 
attributed to Davis, it should be recalled that 

- leal 
i ooo 

whereas out unit of heat is: 

1 JjzL = I cal 

9™ ^80 <lyY^$ 

Thus the two heating rates, of this paper and of Davis, are very nearly 
compatible. The essential results are presented here as meridional 
cross-sections contained in Tables 1 through 6. 

4 

In Table 6b, Davis 1 computations for mean daily Q are based upon 
a climatological atmosphere at rest. It will be seen that Q should 
depend upon the nature of the meridional motion pattern, especially 
the field of vertical motion, in addition to any static distribution 
of radiating material. The term "Balance" is taken from Davis 1 paper 
and refers to balance requirements , on the basis that 

-(fe.R)= q = jjVy • y T+y^*^T-cp-o] 

with the "storage" very small (by several orders of magnitude) both in 
the present study and in Davis 1 . 
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Table 1. Zonal ly-averaged storage rates at lOOmb for period one through nine (1-5 April 1963), in units 
of lCr* gpm sec~l [Conversion factor for ly (gm)~l (day) - !; 1 gpm sec - -*- = 6.2981y (gm) - *- 
(day) - l] 
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Table 2. Zonally-averaged stream-velocity heat transport at lOOmb for periods one through nine 

(1-5 April 1963), in units of 10“^ gpm sec - *- [Conversion factor: 1 gptn sec - '*- = 6.2981y 
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Table 3. Zonal ly-averaged x-~P ot ential advection at lOOmb for periods one through nine (1-5 April 1963) 
in units of 10“^ gpm sec"l [Conversion factor: 1 gpm sec”l = 6.2981y gm”l day"!-] 








1 


r- i 


OsJ | 


r- 


i 


o 


1 


m i 


O' 1 


O' 


i 


o 


1 


— 


1 


O' | 


SO 

O' 




O I 


St | 


00 1 


o 


1 


00 


1 


0- 1 


O' | 


o 


i 


o 


1 


in 


1 


o | 


r— 1 




in i 


Osl | 


r-* | 


ro 


i 


in 


1 


O 1 


ro i 


ro 


i 


r— 


1 


OsJ 


1 


r- i 






1 1 


• | 


• | 


9 


i 


• 


1 


• i 


• i 


• 


i 


9 


1 


• 


I 


• i 


*— < 




55 


1 1 
1 


1 1 
1 




I 

I 




1 

1 


i 

l 


i 

i 


1 


I 

i 


1 


1 

1 




1 

1 


i 

i 


•H 

U 

CL 




1 

I 


1 


1 




I 




1 


i 


i 

« 




l 




1 




1 

a 


i 

a 


< 




1 

1 


1 

1 


1 

I 




i 

i 




1 

1 


i 

i 


1 

i 




l 

i 




1 

1 




1 

I 


1 

i 


m 

1 




1 

1 


1 

r- i 


1 

^ i 


O' 


I 

i 


r** 


1 

1 


i 


l 

r- i 


00 


I 

i 


r— 


1 

1 


ro 


1 

1 


i 

CM | 


»— < i — < 




m i 


ro | 


00 | 


in 


1 


in 


1 


O | 


O 1 


e— 


1 


r- 


1 


O' 


1 


in i 


l 




in i 


r— j 


O | 


St 


i 




1 


OsJ | 


St | 


o 


i 


in 


1 


e— 


1 


CO 1 


0) 




1 1 


• j 


• j 


9 


1 


9 


1 


• | 


• j 


• 


i 


• 


1 


• 


1 


• j 


G 




O 1 
O 1 


1 

r- i 


1 i 

O | 


1 

r- 


l 

i 


ro 


1 

1 


1 1 


1 

r- i 


CM 


i 

i 


l 

CO 


1 

1 


1 

in 


1 

1 


1 1 
0- 1 


•H X) 

G 


o 


1 


CM 1 


O 1 


0- 


i 


ST 


1 


O | 


>o | 


0- 


i 


r- 


1 


•— 


1 


ST 1 


s: i 


CM 


1 


r— | 


OsJ 1 


in 


i 


O 


1 


OsJ | 


O 1 


O' 


i 


o 


1 


*— 


1 


ro i 


£>o B 


1 


1 


• I 


• | 


• 


i 


9 


1 


• 1 


• 1 


• 


i 


9 


1 


• 


1 


• i 


G 60 


IT) 


I 


00 1 


1 


St 


l 


o 


1 


ro i 


in i 


CO 


i 


OsJ 


1 


OsJ 


1 


ro i 


0 

1 i s. 


CM 


1 


CM*— | 


00 | 


O— 


l 


>o*— 


1 


OsJ 1 


roosj i 


O' 


i 


in*— 


1 


'O — 


1 


ro | i 


jG t-h 




o 1 


CO 1 | 


OsJ | 


O' 1 


l 


in 


1 


r— | 


o l 


ro 


i 


«- I 


1 


ro | 


1 


o | 


U 00 




O 1 


CM | 


o | 




1 


OsJ 


1 


00 1 


00 | 


r— 


i 


r— 


1 


r- 


1 


O l 


O' 




1 1 


• 1 


• 


• 


I 


• 


1 


• 1 


• i 


9 


I 


• 


1 


• 


1 


9 | 


<D CM 




m i 


1 1 


CM | 




l 


r— 


1 


1 1 


i 


OsJ 


l 


1 


1 


1 


\ 


CM 1 


C • 




O | 


>o 1 


1 O' 1 


r*- 


i 


U 


1 


St 1 


O 1 


1 o- 


l 


o 


1 




1 


i in i 


O sO 


m 


l 


O I 


0- 1 


OsJ 


l 


o 


1 


St 1 


GO | 


ro 


i 


*— 


1 


O 


1 


O' | 


CO II 


CM 


l 


O 1 


r— | 


r— 


1 


St 


1 


ST 1 


OsJ I 


r- 


i 


o 


1 


r- 


1 


0- 1 


X) 


1 


i 


• 1 


• I 


• 


l 


9 


1 


• 1 


• | 


# 


i 


• 


1 


• 


1 


• | 


o 


o 


i 


O' | 


O 1 


ro 


i 


r— 


1 


OsJ 1 


ST 1 


in 


i 


ST 


1 


st 


1 


O' 1 


•H 1 


ro 


i 


CO | 1 


ST 1 


00 1 


i 




1 


O | | 


00 | 


O' | 


i 


O 


1 


ro i 


1 


CO 1 1 


per 

sec 




m i 


o | 


r— | 


o 


i 


ro 


1 


r*- i 


r- i 


OsJ 


I 


O' 


1 


sT 


1 


00 1 




o r 


O' | 


OsJ | 


r- 


i 


r— 


1 


O' 1 


O 1 


ro 


I 


st 


1 


in 


1 


ST 1 


l. a 




i I 


• 1 


• | 


• 


i 


• 


1 


• i 


• i 


• 


i 


9 


1 


• 


i 


• | 




o I 


OsJ | 


ro i 


*— 


1 


*— 


1 


r— | 


OsJ | 


o 


i 


r— 


1 




i 


ST 1 


O CL 




r- i 


1 O' 1 


1 — 1 


ro 


1 


1 ST 


1 


I CM I 


1 O' 1 


1 ro 


i 


00 


1 




1 


— f'- 1 


U-l 60 


O 


l 


OsJ 1 


o | 


00 


l 


CO 


1 


O' I 


O 1 


OsJ 


l 


00 


1 


0- 


1 


1 h- 1 


X) !— 1 


ro 


I 


in i 


o | 


o 


l 


o 


1 


ro i 


st 1 


o 


l 


CM 


1 


o 


I 


>o 1 


& 


1 


I 


• l 


• j 


9 


l 


9 


1 


• | 


• I 


• 


i 


• 


1 


• 


I 


• 1 


O 


in 


l 


ro i 




r— 


l 


St 


1 


0- 1 


OsJ 1 


r— 


i 


mo 


1 


I 


1 


ro | 


O •• 


ro 


l 




ro | 


St 


I 


in i 


1 


ST | | 


O'*— | 


m | 


l 


ST 


1 


ST 


I 


^J**— | 


*— < S-i 




o | 


OsJ | 




St 


l 


o 


1 


O 1 


o 1 I 




1 


0- 


1 


in 


I 


O' | | 


o 

4-1 




r- i 


00 | 


OsJ | 


OsJ 


1 


o 


1 


OsJ | 


O' I 


OsJ 


i 


in 


1 


r- 


1 


o 1 


cd o 




l i 


• i 


• j 


• 


l 


• 


1 


• | 


• 


• 


I 


• 


1 


• 


1 


• 1 


cd 




in i 


ro | 


1 1 


ro 


l 


CM 


1 




OsJ 1 




l 


1 


1 


ro 


1 


ro | 


4-1 IH 




r- i 


1 o | 


ro | 


ST 


1 


1 — 


1 


l O' 1 


1 o I 


in 


i 


o 


1 


ro 


I 


1 O' | 


u 


in 


i 


O | 


O I 


OsJ 


i 


O 


1 


00 | 


in | 


r- 


i 


CM 


1 


0- 


I 


»— i 


O G 
CL O 


ro 


i 


ro | 


OsJ | 


in 


1 


o 


1 


CM | 


O' 1 


m 


i 


O 


1 


O' 


1 


m i 


CO *H 


1 


i 


• { 


• j 


• 


1 


• 


1 


• j 


• 1 


• 


I 


• 


1 


• 


1 


• i 


G CO 


o 


i 


in i 


ST | 


in 


l 


ro 


1 


1 


0- 1 


m 


i 


CM 


1 


1 


I 


r— | 


Cd Li 


st 


i 


OsJ | 


O 1 1 


o 


i 


m | 


1 


OsJ I 


O I 1 


*— 


i 


CO 1 


1 


in 


1 


CM | | 


Li 0) 

j i s 




in i 


o- 1 


0- 1 


O' 


1 


O' 


1 


o- 1 


00 1 


r- 


l 


o 


1 


o 


1 


O | 


/ % 




h- | 


O | 


in | 


*— 


l 


St 


1 


OsJ | 


st | 


00 


l 


r— 


1 


r- 


1 


00 | 


4-1 O 




l i 


• 1 


• | 


9 


l 


9 


1 


• j 


• | 


• 


l 


• 


1 


• 


l 


• 1 


cd o 




O 1 


OsJ 1 


ST 1 


OsJ 


i 


r— 


1 


1 1 


sT 1 


0- 


i 


f— 


1 




i 


CO 1 


0J . — ■ 




oo 1 


1 ST 1 


r- i 


o 


1 


0- 


1 


>o 1 


OsJ | 


r— 


I 


1 st 


1 


O 


1 


r—r— | 


JG 


o 


1 


ST | 


O' 1 


r— 


l 


O' 


I 


ST | 


St 1 


in 


l 


ro 


1 


ro 


1 


m i 


H 

| 


ST 


i 


OsJ | 


o | 


o 


l 


0- 


1 


O | 


ro | 


ro 


I 


O 


1 


o 


l 


0- | 


cd u 


1 


I 


• | 


• | 


• 


l 


9 


1 


• | 


• j 


• 


i 


• 


1 


• 


l 


• | 


O 0) 


in 


i 


OsJ | 




1 


l 




1 


CM | 


1 




l 




1 


CM 


1 


O 1 


CO 


sT 


l 


00 | 


m i i 


O' 


l 


»— 


1 


St 1 


O' I 


O' 


i 


O 


1 


OsJ 


1 


O | 


4-1 

Li a 
<D rv 




O | 


O 1 


r— | 


in 


I 


St 


1 


r- | 


O' | 


in 


l 


in 


1 


*— 


l 


1 




00 1 


in i 


eg i 


St 


i 


St 


1 


O | 


st | 


ST 


I 


O' 


1 


ro 


1 


O' | 


> 60 




1 1 


• | 


• | 


9 


1 


9 


1 


• j 


• | 


9 


i 


• 


1 


• 


I 


• 1 






in | 


-=3* 1 


O | 


in 


l 


OsJ 


1 


in | 


O 1 


CM 


i 


*— 


1 


»— 


1 


>o | 


X) r-A 




CU | 




O 1 


1 o 


l 




1 


O 1 


1 


•— 


i 


o 


1 


1 CO 


i 


mo i 


<U 1 
60 O 
Cd r— -! 


in 


1 


ro 1 


O' | 


f— 


I 


CM 


I 


00 1 


r— | 


ro 


l 


r*- 


1 


st 


1 


0- 1 




1 


OsJ | 


o | 


o 


I 


O 


I 


O 1 


ST 1 


st 


l 


0- 


1 


*— 


1 


in 1 


Li 


1 


1 


• j 


9 j 


9 


I 


9 


1 


• 1 


• | 


9 


I 


• 


1 


• 


1 


• | 


0) ML 


o 


1 


1 


1 1 


r 


l 


r— 


1 




1 


1 


i 


i 


1 


r— 


1 


1 1 


> o 


m 


1 


m i 


00 1 


00 1 


1 


in 


1 


O 1 I 


O' 1 


0- 


i 




1 


in 


1 


00 | 


cd 

1 CO 

4-1 




m i 


St 1 


O' I 


o 


1 


o 


1 


O 1 


OsJ | 


tn 


I 


o 


1 


in 


1 


r- i 




00 1 


o 1 


O 1 


ro 


l 


in 


1 


00 I 


m i 


CO 


I 


OsJ 


1 


CM 


1 


O | 


i-L *rH 




i i 


• j 


• I 


• 


i 


• 


1 


• 1 


• j 


• 


i 


• 


1 


• 


1 


• j 


rA C 




O 1 


ST 1 


ro l 


OsJ 


i 


>o 


1 


O' | 


r- | 


in 


i 


>o 


1 


CM 


1 


OsJ 1 


cd G 




O' ! 
1 


1 


i 

i 


7 


l 

I 


OsJ 


1 

1 


1 

1 


1 


l 


i 

i 




1 

1 


\ 


I 

1 


in i 
i 


Zor 

in 




J 

/I 


J L 

1 1 




!_ 

l 




JL 

1 


-1— 

! 

1 

1 

J- 

1 

1 

1 




l 




1_ 

1 




1 

T 


i 


H 


/ 1 

/ 1 


1 

1 


l 

1 




l 

l 




1 

1 


1 

1 


1 

1 




1 

l 




1 

1 




1 

1 


l 

GO 1 


<r 


< . 


/ 1 


1 


l 




l 




1 


1 


1 




1 




1 




1 


O-J 1 


0) 




Q 1 


♦— | 


OsJ I 


ro 


1 


St 


1 


in i 


mo | 


I s - 


l 


CO 


1 


O' 


1 


CC<1 1 


T— 1 




CC 1 


1 


1 




1 




1 


l 


1 




1 




1 




1 


CL\- | 






CL | 


1 


1 




I 




1 


l 


1 




1 




1 




1 


1 O | 


cd 




1 


1 


1 




l 




1 


l 


1 




1 




1 




1 


O' h- | 


H 



15 





i 


O | 


ro i 


CO 


1 


O 1 


h- 1 


d 1 


O' | 


in i 


O 1 


ro | 


CO 




O 1 


*" j 


O' 1 


O' 


1 


ro | 


O 1 


h- i 


ro I 


o I 


CM | 


O' l 


\D 




m i 


-d i 


l 


CM 


l 


d 1 


O 1 


m i 


O 1 


O 1 


ro | 


O' l 






I i 


• j 


• i 


• 


1 


• i 


• 1 


• j 


• | 


• i 


• j 


• | 


1—1 




m i 


I i 


i i 




l 


l 


1 


1 


1 1 


i 


1 


1 






m i 


i 


i 




1 


l 


1 


i 


1 


I 


1 


1 


•rl 


* 


i 


i 


i 




l 


l 


1 


i 


1 


i 


1 


l 


u 




i 


l 


i 




1 


1 


1 


i 


1 


i 


j 


l 


a 




i 


I 


i 




» 




1 


i 


1 

i 


1 


1 


l 


< 




1 

i 


CM 1 


i 

CO | 


CO 


1 

1 


I 

d | 


1 

GO | 


l 

00 1 


l 

00 1 


i 

•O 1 


1 

ro i 


l 

ro i 


m ■ — i 
1 i 




in i 


O | 


p 1 


CM 


1 


CO 1 


O I 


CO | 


CM 1 


O 1 


O 1 


ro I 


1 1 




in i 


in | 


r— j 


ro 


1 




ro | 


p j 


m i 


CM 1 


-d 1 


O' | 


^ B 




I i 


• j 


• I 


• 


I 


• j 


• j 


• j 


• i 


• | 


• | 


• j 


bO 




O 1 


1 


1 1 


1 


1 


1 


1 




1 


1 


1 


CM 1 


a) 

1 




O 1 


C7> | 


O 1 


o 


1 


O 1 


CM 1 


d 1 


G | 


O' | 


00 | 


O 1 


•H I 


o 


i 


p- | 


O' 1 


o 


l 


o | 


p j 


o 1 


h- 1 


CM | 


d- i 


in | 


a >> 


CVJ 


i 


CM 1 


r— | 


in 


1 


*0 1 


CM 1 


O 1 


O' | 


| 


o | 


ro i 




1 


i 


• | 


• 1 


• 


l 


• 1 


• | 


• 1 


• I 


• j 


• j 


• | 


XI X) 


UO 


i 


CO 1 


1 


d 


1 


O 1 


ro 1 


in i 


CO | 


CM | 


CM | 


ro | 


60 

d >% 


CM 


i 


o*— I 


p | 


o— 


1 


cr*p* i 


CM l 


r- CM | 


CM | 


CM*— | 


Op- | 


ro i | 


O 




O 1 


PO | | 


CM | 


d 1 


1 


O 1 


o l 


m 1 


m i 


CM | | 


CO 1 | 


1 


M 00 




O 1 


ro | 


r- i 


o 


l 


ro 1 


p | 


O I 


O 1 


CM | 


p j 


1 


X2 ON 




1 1 


• j 


• i 


• 


l 


• 1 


• j 


• I 


• i 


• | 


• | 


• i 


U CM 




in i 


1 1 


p | 




1 


p | 




1 1 


ro | 


p j 




O 1 


• 

n \ in 




"O 1 


m i 


1 g | 


CM 


1 


O' | 


i r- 1 


O' | 


1 d i 


1 d- | 


1 ro i 


1 O 1 


\U MJ 

rj 


to 


i 


ro i 


O' 1 


CO 


1 


CO 1 


O' | 


GO | 


O | 


m i 


O' | 


-d l 


o II 


CVJ 


i 


in I 


p I 


p- 


1 


d 1 


ro | 


CM | 


O I 


G | 


o 1 


d* 1 




1 


I 


• | 


• 1 


• 


1 


• 1 


• j 


• | 


• 1 


• | 


• j 


• | 


XJ 1-1 


o 


i 


O' | 


o 1 


ro 


1 




CM | 


d i 


in i 


-d | 


-d i 


O' I 


O 1 


ro 


I 


O 1 1 


O | 


d- 1 


1 


d | 


*- 1 | 


in | 


— i i 


o | 


ro 1 1 


*- 1 1 


•rH KJ 

J-l <1) 




in i 


r— 1 


ro | 


h- 


1 


O | 




d | 


in | 


ro | 


O | 


in i 


<D CO 




O 1 


ro l 


O' | 


in 


1 


O 1 


d i 


00 | 


p* j 


O | 


r— j 


O' | 


a 




I l 


• i 


• j 


• 


1 


• j 


• j 


• 1 


• i 


• j 


• j 


• j 


S 




O 1 


ro | 


ro i 


p- 


1 




CM | 


CM | 


O | 


CM 1 


p | 


-d | 


u a. 




h- i 


1 CO 1 


1 CO l 


d 


1 


I d i 


i in i 


1 in i 


1 ro i 


O 1 


O | 


p-ro | 


O 00 

[ ) j 


o 


l 


p l 


-d | 


ro 


l 


r- i 


ro i 


O 1 


in i 


in i 


h- i 


1 in | 


H 


ro 


l 


O | 


O 1 




l 


m i 


CM | 


ro | 


o I 


CM | 


^ l 


ro i 


X) 


1 


l 


• 1 


• 1 


• 


1 


• j 


• | 


• i 


• i 


• 1 


• j 


• j 


s 


m 


l 


ro | 


p- 1 


p- 


1 


d 1 


h- 1 


CM 1 


r— j 


O 1 


1 l 


ro | 


G •• 


ro 


i 


O | 


JT | 


h- 


l 


O | 1 


d i | 


in*— i 


r- i l 


d- | 


G 1 


f^-»— | 


O M 

T— 1 O 




O 1 


r— j 


O | 


o 


1 


ro i 


ro I 


d | | 


m i 


CM 1 


ro i 


O | I 


u 




r- i 


CM | 


p- | 


p- 


1 


d i 


d I 


ro i 


O 1 


CM 1 


ro l 


ro i 


4J o 




i i 


• 1 


• j 


• 


1 


• I 


• j 


• l 


• | 


• | 


• | 


• | 


CCJ Cti 




in i 


-=r i 


1 | 


ro 


1 


CM 1 




ro | 


I i 




ro | 


O | 


•IM 




h- i 


1 CM 1 


ro | 


in 


1 


1 in I 


1 ro | 


1 *o | 


-d | 


1 CM | 


ro | 


1 m i 


CO 

<D d 


in 


i 


-d | 


ro I 


in 


1 


ro | 


CM | 


O' | 


r- i 


O | 


d - | 


p i 


u o 


m 


1 


ro | 


CM 1 


in 


1 


O | 


ro | 


CO 1 


m i 


O | 


O' | 


ro | 


cij »H 


\ 


1 


• i 


• | 


• 


1 


• j 


• j 


• j 


• j 


• j 


• j 


• j 


U CO 


o 


l 


in i 


-d | 


in 


1 


ro | 


1 


r- i 


in i 


CM | 


1 1 




l-i 


d* 


i 


CO | 


r- 1 l 


d 


1 


ro | | 


00 1 


m i i 


CM | 


ro | | 


ro 1 


ro | | 


00 CJ 
d > 




in i 


o | 


CJ | 


O' 


1 


ro i 


CO 1 


in | 


O 1 


O 1 


i 


O | 


•h d 




r- i 


f— 1 


d" | 


o 


1 


CM | 


d i 


ro l 


CO | 


CM | 


ro I 


-d | 


4J O 




I i 


• 1 


• | 


• 


1 


• | 


• i 


• | 


• i 


• | 


• | 


• | 


d CJ 




O 1 


CM 1 


d \ 


CM 


1 




l 


d 1 


r- i 


p | 


r— | 


r- i 


CJ 1 — • 




CO | 


1 m i 




h- 


1 


h- 1 


i 


ro 1 


ro 1 


I ro | 


p- j 


*-cr i 


i-G' 


O 


1 


>o 1 


I s - | 


CM 


1 


r— | 


O' | 


O' | 


CM | 


G | 


d- | 


CM | 


o | 


-d 


1 


p 1 




p- 


1 


O | 


O' | 


CM 1 


ro | 


o | 


G | 


r— | 


•H O 


i 


1 


• i 


• j 


• 


1 


• | 


• 1 


• l 


• | 


• 1 


• | 


• 1 


4J CJ 


in 


1 


CM 1 


p- | 


1 


1 


1 


p- 1 


1 


1 


1 1 


CM | 


O 1 


d co 


d* 


1 


O 1 


d 1 1 


ro 


1 


O 1 


CM | 


o 1 


d- | 


O 1 


CM | 


CM | 


n 

d g 




O 1 


CO 1 


d 1 


o 


1 


CO | 


O | 


CO 1 


CO | 




ro | 


O | ' 


•H (X 




CO 1 


-d 1 


p | 


ro 


l 


r- i 


O | 


d | 


rr) i 


CM 1 


O | 


r- i 


XJ 00 




1 1 


• | 


• | 


• 


l 


• i 


• j 


. • I 


• j 


• j 


• j 


• i 






in i 


.d 1 


O I 


in 


1 


CM 1 


m i 


O I 


CM | 


CM j 


p | 


r- i 


XJ 1—1 

1 




CO | 


O' I 


X-> I 


i 


i 


«— IM | 


O' 1 


*— LH | 


•— i 


G 1 


1 n> 1 


r*)*— i 


1 

00 o 


n 


1 


m i 


p | 


ro 


l 


in i 


O 1 


h- 1 


p i 


ro | 


in | 


CO 1 


cd i—i 


d- 


1 


in i 


O | 


in 


1 


r- i 


O 1 


d | 


in i 


r- i 


*— j 


O 1 


u 


I 


1 


• j 


• j 


• 


1 


• j 


• i 


• j 


• i 


• i 


• i 


• 1 


CJ 14-1 


o 


1 


l 


1 1 


p- 


1 


p | 


p- | 


1 


l l 


i i 


r— 1 


1 


> o 

rH 


in 


1 


-d i 


ro I 


d 1 


1 


r- i 


O 1 1 


r- 1 


O 1 


ro 1 


CM 1 


00 I 


1 co 




in | 


o | 


-d 1 


d* 


1 


CO 1 


CM | 


C0 1 


O 1 


p | 


CO | 


CO 1 


>> 4J 




CO | 


O 1 


CO 1 


O 


1 


O 1 


in i 


00 | 


00 1 


-d | 


O | 


CM 1 


rH »H 




i i 


• I 


• I 


• 


1 


• i 


• i 


• 1 


• 1 


• j 


• I 


t | 


i— 1 1 d 




O 1 


in i 


ro 1 


CM 


1 


r- i 


O 1 




in | 


O 1 




in i 


cd d 




O' 1 


p i 


1 


p- 


1 


CM 1 


p 1 


*— 1 


l | 


1 


1 1 


in | 


o d 




i 


i 


l 


1 


l 


i 


1 


\ 


1 


1 


1 


1 






i 


i 


J 




1 


i 




A. 


J 


J. 


4 


l 








i 


1 




l 


T 


__ .. J- 


i 


1 


1 


1 


1 




H 


/ 1 
/ i 


l 

1 


1 

1 




1 

1 


1 

1 


1 

1 


i 

i 


1 

1 


i 

i 


1 

1 


l 

00 1 


m 


< 


/ i 


1 


1 




1 


1 


1 


i 


1 


i 


1 


a_j i 


CJ 


-jy 


O 1 




CM 1 


ro 


1 


d 1 


in i 


O | 


h- 1 


co i 


O' 1 


C£<L 1 


i— i 




CC | 


1 


I 




1 


1 


i 


1 


1 


I 


1 


CLh- 1 


JO 




a. i 


1 


1 




1 


1 


1 


1 


1 


I 


1 


1 G 1 


cd 




I 


l 


1 




1 


1 


i 


1 


1 


I 


1 


O'H* 1 


t-* 



16 



o 

r 

ur\ 

03 


Sr 


s ! 


i £3 


oj 

cG 

VO 


s 


LQl 

u 0 

LTV 


UN 

? 

CO 


£ 


c r 
o 


« 


cO 

<r 


f 




9 > 

j 

if* 

r* 


oo 

CM 


t£> 

l' 


cO 

C 4 


cr 

oi 


cM 

i 




in 

rr 

o 

h 


Sj 


CT- 

<M 

r 


O" 

c 4 


$ 


IS 

r 


VO 

o 

# 


o 

«r 

iS 


oO 

o 

r 


ri 

CO 

r 


o 


r 


ES 


r- 

\ 


$ 

s 


■>£ 

r 


r 


O 

* 


i ' 


rrj 


f- 


0 

1 

u* 

UN 


i « 

\ 


•<0 

O 


s 
1 ' 


o-l 

r 

• 


n 

1 r 


O 

IA 

8 


cC>i 

o 

1 - 


K} 


; oo 

o 

H 


lO 

o 

# 


o 

7 


o 
1 # 


o 

* 


C*i 

o 


in 

o 


v 9 

O 

r 




td to 


\p 

t 

o 


a 


0 

1 • 


N- 

G 

< ' 




Si 
1 1 


r- 

o 


0 

1 

lA 

K* 


•3 

o 


o 

o 

r 


O 

t* 


j 


rs! 

X 


1 ' 


lA 

•0 

$ 


to 

o 


«N 

O 


o 

r 


9 

1 


*0 

<r 


r 


S 

A 


O 

o 

+' 


i. ^ 

■ o 


o 

o 

l ' 


g 

T 


cr 

cr 


c- 

cr 

r 


VI 

I 

s 


Q 

V 


ij o 


o 
0 
1 V 


vA 

? 


U? 


in 

U 3 

i ’ 




N 

*C 

* 

4 

Sb 

J 


i 

- 9 

i > 

4 ^ 


■ *- 
► 


b 

^ 3 
d V 


§ 

Ui 

o 

^-7 

‘1 


it 

Jir 

J*Q 

Cr 



0 

r 

u> 

o 


sP 

fC 


to 

« 


i 0 

cvl 


X 

cJ 

« 


S3 


cr 

?- 


in 

<5 

it 


o 


X 

# 


S 


S 3 

rM 


oO 

cy 


t * 


o 

* 

U* 

U 5 


l/l 

• 


in 


>”> 

f 4 


tsl 

Cvi 


O 

c^i 


r- 


lA 

« 

O 

to 


cr 

o 


- 


CM 

cM 




cr 

# 


l' 


0 
na 

1 

rr 


$ 


v£ 


Q 

Cvl> 

* 


cr 


VJP 


in 
t • 


tA 

Ck 


o 

o 

JL 




cO 




to 

« 


± 

r 


0 
<r 

1 

iA 

to 


0 


c£ 


•£ 


u} 


# 


r 




i* 


0 


X 


K} 


o 

# 


<T“ 

o 

r 


s 

a 


O 


& 


c 4 

0 


o 


CO 

o 


r* 

o 

r 


t £5 

o 


rO 

O 


£ 

$ 


O 


r~ 

o 


C 


^r 

o 

r 




* £ 
4 

sy* 


r° 

§c 

• 0 

✓> 


i— 

3 ‘o 
2 v 
r (>■ 

s° 

C 

WN 


vU O 

a o 

wi 

20 

H 


OJ 

s 

<t 

c 

cO 


ii 

3 ! 

‘Ti 

cr 



u 

u 

o 

a 

cn 

a 

oj 

u 

4-> ^ 



o3 

03 



oj 

cs 

o 

•H 

xj 

U 

03 

B 



u 

u 

4) X> 

> 4) 

00 
03 
U 

03 < 
> 
a 3 



C 

0 

*(— I 

u 

a 

a; 

> 

x) 

03 

1 

X 



- rC 

c 

o u 



•* in 
00 -H 

a > 

•H 03 
4 -> Q 
o3 

03 4-1 



4_» 

0 
0 ) 
> 
X) 
o3 

1 



a; 

cn 

o 

-u 



o3 
£i 
03 
•H X 
XJ •*-» 

•iH 

$ 

r—l 

•H IX 
03 XI 



03 XJ 

00 

03 03 

U rC 

O 4 -> 

4-> 

cn o 

u 



3 

u 

in 



rC 

u 



O cn 4-i 
C O 
c n o 
4) »H *0* 
00 
03 
U 
4) 

> 

03 

3-i 
P 

0 
X 

1 

<r 



3 4-1 
& O 

•rl 

J-« c 
u o 
C cn 

O 'H 
O. 3-1 
03 

n 9- 
.*h E 

4) O 
CNI’ & O ‘ 



x> /-N 

a 

o3 ^ 



vO 

03 

i— l 

03 

H 



17 



5. 



Conclusions: stratospheric circulations and associated physical 



mechanisms 

As already noted, the zonally-averaged storage term is very small. 
In addition the stream- and potential-advections , averaged over five- 
degree latitude bands are smaller by approximately two orders of magni- 
tude than [Q] . Thus to a good approximation, valid everywhere at 
least as to sign, 



’q]= --c P [riG] 



and is of the proper order of magnitude as well. 
By the equation: 



(17) 



it follows that [^qq] >0 occurs when [^. 'V] represents convergence. 
Now convergence occurs near the pole , since the expression for the 
zonally-averaged divergence 




(18) 



with the second term (which represents the convergence of meridians) 
becoming dominant [Martin, 4] at latitudes (|) > 75N. At lower latitudes 

[ >' = (V^ooU = [3%^ 



is a valid approximation; that is, the spherical divergence is only 
slightly affected by convergence of meridians. Hence, where (^iQgj i- s 
negative, there is divergence^ if^Q^jis P osit ive, there is convergence* 
Consequently at lOOmb, we have centers of convergence at the pole and 
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latitude 35, and of divergence at latitude 65. 

In inferring a stratospheric circulation, one may start with an 
interhemispheric meridional circulation proposed by Kellogg [6], 
especially applicable, in some of its aspects, to the winter-summer 
reversal in the stratosphere. The octagonal grid, together with the 
associated boundary conditions, is not able to resolve the effect of 
the cross-equatorial flow to at least latitude 20N latitude. However, 
because of the decrease of mass in the winter stratosphere, cross- 
equatorial flow must occur near the time of the vernal equinox. This 
evidently is responsible for the convergence center near latitude 35N. 

Likewise with the return of the sun to the Northern Hemisphere, 
ozone regeneration goes on rapidly, followed by ozone-heating, 
especially at 45 km. According to Pressman [7]-, an ozone radiative 
heating maximum proceeds progressively from 30N latitude to the pole 
between the dates 15 March and 15 May. For time of this study, 1-5 
April, 1963, one might expect maximum ozone heating at 45 km at or near 
60N latitude. This effect is so marked as to cause a strong forced 
outdraf t (V‘V>0 ) at 60N latitude, and at an elevation of 45 km. In 
order that three-dimensional mass divergence be within reasonable 
bounds, there must be a sympathetic updraft from below. This is the [co] 
measured at the lOOmb level and presumably exists at levels between 
lOOmb and 45km. Also according to Martin’s model, sizeable co-values of 
the same sign still persist at the 200mb level, but are smaller than 
co^QO units of cm sec”l. The level of "return flow" in these cells 
cannot be designated from presently known studies. 
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fvC-U VA 11-3 




The ’’maximum" convergence in the vertical and associated downdrafts at 
the pole at this particular time of ozone activity should then be 
viewed as a "spill-over” effect from 65N latitude. The maximum con- 
vergence at or near latitude 35 North, which mirrors a zonally-averaged 
100-mb maximum values of a) at this same latitude is postulated as the 
primary result of cross-equatorial flow. Considering the level at which 
the greatest polar-vortex intensification in the Southern Hemisphere 
occurs, the level of the 35N latitude maximum of convergence has also 
been placed at 45km. 

Fig. 4 together with the above discussion "explains" the mean 
meridional flow pattern in the stratosphere during the period of this 
study. One may offer as conjecture that as the stratospheric summer 
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advances (say by June 15), the maximum convergence will have advanced 
to 45N latitude while the zone of maximum divergence associated with 
most effective ozone-heating near 45km advances to the pole, there to 
become associated with divergence at higher elevations. 



moderate cooling at every latitude due to predominance of long-wave 
radiation over heat sources. Two features by which April, 1963 
appears to differ from climatology are shown in Figs. 5 and 6, both 
taken from analyses of the Free University of Berlin [8]. The latter 
figure shows a pre-existing zone of vortices at latitude 60 North, 
as of 1 April, 1963, 1200 GMT. The other figure mentioned above shows 
that, as compared with a completely different April (1962), pressure 
was higher over the pole but lower than in 1962 in a vortex ring at 
60N latitude. Fig. 6 shows a possible reason for this: strong baro- 

clinicity on the equatorward side of the vortices, with f zonal 
generation' _ | / / 



computed over 45-70 degrees North latitude favoring the continuation of 
such a baroclinic zone and westerly winds. In fact, according to the 
analyses of the Free University of Berlin, an easterly regime with an 
associated polar anticyclone did not occur until May 15, 1963. 

Finally, how may one reconcile the zonally-averaged equation: 



so that both horizontal flows and storage are relatively unimportant 



As contrasted with the Q of this study, Davis' results show only 



G 



z 





(19) 
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Fig* 5 . lOOmb height changes, April 1963 - April 1962 from analyses of the 
Free University of Berlin 
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Fig. 6. lOOmb analyses, OOGMT , 1 April 1963 from analyses'- of the Free University 
of Berlin 
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as determined here. Equation (19) implies that with >0 (down 

draft motions) , £qJ<0 is the result. Such a result is reasonable on 
physical grounds by noting that descending air is warmed dry adiabati- 
cally, and in addition one may expect some increased diabatic solar 
heating acting upon an ozone increase in lower stratospheric elevations 
However, one immediate opposing result is that the CO 2 in the warmed 
air will emit at a greater rate, tending to re-establish the heat 
balance. The same argument applies in reverse for updraft motions. 
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